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We review the status of Monte Carlo generators presently used for simulations of the large-angle Bhabha 
process at electron-positron colliders of moderately high energy (flavour factories), operating at centre-of-mass 
energies between about 1 GeV and 10 GeV. It is shown how the theoretical accuracy reached by present Bhabha 
programs for physics at flavour factories is at the level of 0.1% and, therefore, comparable with that reached about 
a decade ago for luminosity monitoring through small-angle Bhabha scattering at LEP. 



1. INTRODUCTION 

At modern electron-positron (e + e~) accelera- 
tors built for a deeper and deeper understanding 
of the fundamental interactions between elemen- 
tary particles, the machine (integrated) luminos- 
ity can be derived precisely by means of the re- 
lation J Cdt = iVobs/oth, where iV bs and <7th are 
the number of events and the theoretical cross 
section of a given reference process, respectively. 
Because of the latter relation, the reference pro- 
cess must be a reaction with clean topology, high 
statistics and calculable with high theoretical ac- 
curacy, to maintain small the total luminosity er- 
ror given by the sum in quadrature of the relative 
experimental and theoretical uncertainty. 

For example, at high-energy accelerators 
LEP/SLC running in the '90s around the Z pole 
to perform precision tests of the Standard Model, 
the process of e + e~ production (Bhabha scatter- 
ing p]), with the final-state leptons detected at 
small scattering angles, was used because dom- 
inated by the electromagnetic interaction and, 
therefore, calculable in perturbation theory, at 
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least in principle, with very high theoretical ac- 
curacy. A total (experimental plus theoretical) 
precision of ~ 0.05 -j- 0.1% was achieved at the 
end of LEP/SLC operation [2"I3"H] . thanks to the 
work of different theoretical groups and the ex- 
cellent performances of precision luminometers. 

At presently running e + e~~ colliders of mod- 
erately high energy (between about 1 GeV 
and 10 GeV), such as the ^-factories DA$NE 
and VEPP-2M, the charm-factories CESR and 
BEPC, as well as the B-factories KEK-B and 
PEP-II (globally denoted as flavour factories), 
the normalization reaction primarily used is the 
large-angle Bhabha process [5]. Actually, at all 
flavour factories, the final-state leptons are de- 
tected at wide scattering angles, within a typi- 
cal central acceptance region of ~ 30° 4- ~ 150°, 
because of the absence of dedicated luminosity 
counters, for example, at small scattering angles. 
A precision luminosity measurement at flavour 
factories is of utmost importance to perform ac- 
curate measurements of the e + e~ — > hadrons 
cross section, which is, in turn, a key ingredient 
in high-precision calculations of the running of 
ciqed and lepton g — 2. 

It is worth noting that high theoretical accu- 
racy and, in particular, comparison with pre- 
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cise data require the development of sophisticated 
Monte Carlo (MC) programs, including radiative 
corrections with a relative precision at the per 
mille level. In this review, only the large-angle 
Bhabha generators used at present e + e~ collid- 
ers will be considered. A description of the small- 
angle Bhabha programs available in the literature 
and used at LEP can be found in [6]. 

2. LARGE ANGLE BHABHA GENERA- 
TORS 

Since pure weak corrections are completely ir- 
relevant in the energy range explored by flavour 
factories, the dominant part of radiative correc- 
tions, i.e. leading-log QED contributions, can be 
kept under control by means of different theoret- 
ical methods, such as QED Structure Functions 
(SF) [7] and Yennie-Frautschii-Suura (YFS) ex- 
ponentiation [H], as done in many applications to 
precision physics at LEP [3] . For instance, in the 
generator BabaYaga v3.5 |9ll0j , a pure QED Par- 
ton Shower (PS) is implemented, through a MC 
solution of the QED DGLAP equation satisfied 
by the electron SF in the non-singlet approxima- 
tion. Once the electron SF is numerically derived, 
the QED corrected cross section can be simply ob- 
tained by dressing the hard-scattering cross sec- 
tion of the process under consideration with a SF 
for any charged leg. A remarkable advantage of 
this universal approach is that the four-momenta 
of the final-states products (charged particles and 
photons) can be exclusively generated through 
the shower cascade, as in the case of the QCD 
PS programs for the momenta of quarks and glu- 
ons. 

A comparison between BabaYaga predictions 
and data published by the KLOE Collaboration 
for the measured energy distribution of the e + e~ 
pairs produced in the large-angle Bhabha process 
at DA$NE is reported, from [IT], in Fig. [TJ show- 
ing a nice agreement between data and theory. In 
spite of that, it is clear that an obvious drawback 
of a pure PS generator like BabaYaga v3.5 is that 
it lacks the effect of 0(a) non-logarithmic con- 
tributions and, therefore, the corresponding theo- 
retical accuracy is limited to that of a leading-log 
(LL) approximation. However, the physical pre- 




Figure 1. KLOE large-angle Bhabha data for 
the energy distribution of the e + e~ final state. 
From [IT] . 



cision of the formulation underlying such a pro- 
gram can be substantially improved from about 
a per cent accuracy down to the per mille level 
by performing a matching of the exact next-to- 
leading (NLO) corrections with the multiple pho- 
ton effects taken into account by the PS. This 
matching procedure, which is described in detail 
in [12], allows to obtain an O(o) corrected cross 
section that coincides, by construction, with the 
exact NLO result, avoiding double counting and 
preserving resummation of leading contributions. 
This theoretical formulation is on the grounds of 
an improved version of BabaYaga v3.5, denoted 
as BabaYaga@NLO, having a precision tag of 
0.1%. Needless to say, thanks to the matching al- 
gorithm, the theoretical error of BabaYaga@NLO 
is shifted to 0(a 2 ) contributions of the next -to- 
next -to-leading (NNLO) perturbative expansion. 

Theoretical approaches rather similar, at 
least in their basic ingredients, to that of 
BabaYaga@NLO are implemented in the MC pro- 
grams MCGPJ [13] and BHWIDE [IT], both 
used, as BabaYaga@NLO, at flavour factories for 
simulations of the Bhabha process. The former is 
a generator realized by a Dubna-Novosibirsk col- 
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Table 1 

Relative size of different sources of correction (in per cent) to the large-angle Bhabha scattering cross 
section for realistic luminosity set up at <£>- and i3-factories. From [T2] . 



set up 


a. 


b. 


c. 


d. 


5 a 


-13.06 


-17.16 


-19.10 


-24.35 


rnon-log 


-0.39 


-0.66 


-0.41 


-0.70 


<5ho 


0.43 


0.93 


0.87 


1.76 


$a 2 L 


0.04 


0.09 


0.06 


0.11 


<5vp 


1.73 


2.43 


4.59 


6.03 



laboration and includes exact 0(a) corrections 
supplemented with higher-order LL effects taken 
into account through collincar SF, resulting in 
a theoretical precision estimated by the authors 
to be better than 0.2%. The latter is a MC 
code developed by the Krakow-Knoxville group 
at the time of LEP operation and implements ex- 
act 0(a) corrections matched with the resumma- 
tion of soft and collinear logarithms via YFS ex- 
ponentiation. According to [13], the precision of 
BHWIDE is conservatively estimated about 0.5% 
for LEP1. However, since the theoretical ingredi- 
ents of BHWIDE are very similar to the formu- 
lation of both BabaYaga@NLO and MCGPJ, it 
is reasonable to assume that its theoretical accu- 
racy for physics simulations at flavour factories is 
at the level of 0.1%, as discussed in the following. 

Concerning the theoretical precision of all the 
above generators, it is important to empha- 
size that the bulk of the most important sub- 
leading 0(a 2 ) corrections, namely a 2 L pho- 
tonic contributions enhanced by infrared log- 
arithms, where L = ln(Q 2 /m 2 ) is the large 
collinear logarithm, is effectively incorporated 
in the tools BabaYaga@NLO, MCGPJ and BH- 
WIDE by means of factorization of 0(a) non-log 
terms with the leading 0(a) contributions taken 
into account in the PS, collinear SF and YFS ex- 
ponentiation approaches, as demonstrated in |15j . 

3. NUMERICAL RESULTS 

To get an idea of which corrections are rele- 
vant to achieve a per mille precision in simulations 
of the Bhabha scattering at flavour factories, we 
show in Tab.Q]the relative effect of various contri- 



butions to the large-angle Bhabha cross section, 
when considering typical selection criteria at <fr— 
(set up a. and b.) and i3-factories (set up c. and 
d.) (see [12] for details). 

From Tab. [TJ it can be seen that 0(a) cor- 
rections decrease the Bhabha cross section of 
about 15% at the ^-factories and of about 20- 
25% at the B-factories. Within the full set of 
0(a) corrections, non-log terms are of the or- 
der of 0.5%, almost independently of the centre- 
of-mass (cm.) energy, as expected, and with a 
mild dependence on the angular acceptance cuts, 
as due to box/interference contributions. The 
effect of higher-order corrections due to multi- 
ple photon emission is about 0.5-1% at the <&- 
factories and reaches 1-2% at the i?-factories. 
The contribution of (approximate) 0(a 2 L) cor- 
rections is not exceeding the 0.1% level, while the 
vacuum polarization increases the cross section 
of about 2% around 1 GeV and of about 5-6% 
around 10 GeV. Concerning the latter correction, 
the non-perturbative hadronic contribution to 
the running of a is included in BabaYaga@NLO 
both in the lowest-order and one-loop diagrams 
through the HADR5N routine [16], that returns 
a data driven error, thus affecting the accuracy 
of the theoretical calculation. Analogous results 
about the size of radiative corrections have been 
obtained recently [17] for the process e + e~ — > 77, 
also of interest for precision luminosity studies at 
flavour factories. As a whole, these results indi- 
cate that both exact 0(a) and higher-order cor- 
rections (including vacuum polarization) are nec- 
essary for 0.1% theoretical precision. 
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Figure 2. Tuned comparisons between the 
results of BHWIDE, BabaYagaONLO and 
BabaYaga v3.5 for the electron energy distribu- 
tion in the large-angle Bhabha process at PEP- 
II. Courtesy of A. Denig and A. Hafher. 



4. TECHNICAL AND THEORETICAL 
ACCURACY 

4.1. Technical precision: tuned compar- 
isons 

A typical procedure followed in the literature 
for establishing the technical precision of MC 
generators is to perform tuned comparisons be- 
tween independent predictions, using the same 
set of experimental cuts. An example of such 
tuned comparisons is given in [12] . where it is 
shown that the agreement between the predic- 
tions of BabaYaga@NLO and BHWIDE at the 
^-factories is well below 0.1%, and that also the 
agreement between BabaYaga@NLO and LAB- 
SPV, which is a benchmark code by our group 
with a formulation based on collinear SF very 
similar to MCGPJ, is very good, below the 0.1% 
level. This level of agreement, together with 
further considerations about the size of two- 
loop corrections discussed in the next subsection, 



is the reason why in the latest publication by 
KLOE Collaboration about the measurement of 
the hadronic cross section at DA<I>NE [TH] the rel- 
ative uncertainty assigned to theory in the lumi- 
nosity measurement is now 0.1%, resulting in a 
total luminosity error of 0.3%. 

Similar comparisons have been performed be- 
tween the results of BabaYaga@NLO and BH- 
WIDE by A. Denig and A. Hafner of BABAR Col- 
laboration, in the presence of realistic selection 
cuts for luminosity at PEP-II. Their studies show 
that the two generators agree at the 0.1% level, 
both for integrated and differential cross sections, 
in the physical and statistical significant regions. 
An example of such a comparison, showing the 
predictions of BHWIDE, BabaYaga@NLO and 
BabaYaga v3.5 for the electron energy distribu- 
tion (upper panels), as well as the relative differ- 
ences between the results of BHWIDE and the 
two BabaYaga versions (lower panels), is given in 

Fig. m 

4.2. Theoretical precision: comparisons 
with two loop calculations 

In order to assess the physical precision of the 
generators, the methods typically used are i) to 
compare with 0(a 2 ) calculations, if the latter - 
as in the case of Bhabha scattering - are avail- 
able in the literature [19|20|21|22|23|24j ii) to es- 
timate the size of unaccounted higher-order con- 
tributions. 

Concerning point i) and considering, for def- 
initeness, the generator BabaYaga@NLO, the 
strategy consists in deriving from the general for- 
mulation the 0(a 2 ) cross section, which can be 
cast in the following form 



'sv 



°SV.H + °HH 



(1) 



where, in principle, each of the above 0(a 2 ) con- 
tributions is affected by an uncertainty, to be 
properly estimated. In Eq. ([1]), the first contri- 
bution is the cross section including 0(a 2 ) soft 
plus virtual corrections, whose uncertainty can 
be evaluated by comparison with the available 
NNLO calculations. The crgy of BabaYaga@NLO 
has been compared, in particular, with the cal- 
culation by Penin [21] , who computed the com- 
plete set of two-loop virtual photonic corrections 
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Figure 3. 

BabaYaga@NLO and the 0{a 2 ) calculations 
(photonic) and [55] (Np = f), as a function of 
the infrared regulator. From [T2] . 



Figure 4. The same as Fig. 3, as a function of a 
fictitious electron mass. From H2l. 



in the limit Q 2 3> m 2 supplemented by real soft- 
photon radiation up to non-logarithmic accuracy, 
and the calculations by Bonciani et al. [22] , who 
computed two-loop fermionic corrections (in the 
one-family approximation) with finite mass terms 
and the addition of soft bremsstrahlung and real 
pair contributions. 

The results of such comparisons are shown in 
Fig. [3] and in Fig. 2] for set up a. at the <&- 
factories. In Fig. [3] 8a is the difference between 
tJgy of BabaYaga@NLO and the cross sections 
of the two Oia 2 ) calculations, denoted as pho- 
tonic (Pcnin) and Np = 1 (Bonciani et al.), as 
a function of the logarithm of the infrared reg- 
ulator e. It can be seen that the differences are 
given by flat functions, demonstrating that such 
differences are infrared-safe, as expected, as a 
consequence of the universality and factorization 
properties of the infrared divergences. In Fig. [4j 
8a is shown as a function of the logarithm of a 
fictitious electron mass and for a fixed value of 
e = 10~ 5 . Since the difference with the calcu- 
lation by Penin is given by a straight line, this 
indicates that the two-loop soft plus virtual pho- 
tonic corrections missing in BabaYaga@NLO are 
0(a 2 L) not infrared-enhanced contributions. On 
the other hand, the difference with the calcula- 
tion by Bonciani et al. shows that the fermionic 



two-loop effects missing in BabaYaga@NLO are 
dominated by 0(a 2 L 2 ) contributions. It is im- 
portant to emphasize that, as shown in detail 
in [T2], the sum of the differences with the two 
0(a 2 ) calculations does not exceed the 1.5 x 10 -4 
level, for all the considered set up at <&- and B— 
factories. The second term in Eq. fl} is the cross 
section containing the one-loop corrections to sin- 
gle hard bremsstrahlung and its uncertainty can 
be estimated by relying on partial results exist- 
ing in the literature. Actually, the exact pertur- 
bative expression of Cgy H is not available yet for 
full s + t Bhabha scattering, but, using the re- 
sults valid for small-angle Bhabha scattering 25J 
and large-angle s-channel processes [5^, the rel- 
ative uncertainty of BabaYaga@NLO in the cal- 
culation of (Tg V H can be safely estimated at the 
level of 0.05%. The third contribution in Eq. (JTJ) 
is the double hard bremsstrahlung cross section, 
whose uncertainty can be evaluated by compar- 
ison with the exact e + e~ — > e + e _ 77 cross sec- 
tion. As shown in [12], the differences registered 

2 

between a^ R as in BabaYaga@NLO and the ex- 
act calculation are really negligible, at the 10 -5 
level. 

Summing all the results for the various sources 
of uncertainty, it turns out that the total theoret- 
ical error in BabaYaga@NLO is ~ 0.1% , when 
also including the uncertainty due to the running 
of a as returned by the HADR5N routine and the 



6 



G. Montagna 



contribution, at a few 0.01% level, of light pairs 
radiation [Hj , still missing in BabaYaga@NLO. 

5. CONCLUSIONS 

During the last few years, there has been a sig- 
nificant progress in reducing the theoretical un- 
certainty in Bhabha generators used at presently 
running e + e~ colliders down to 0.1%. Exact 
0(ot) and multiple photon corrections are nec- 
essary ingredients to achieve such a precision. 
These corrections are implemented in three gener- 
ators (BabaYaga@NLO, BHWIDE and MCGPJ) 
for the large-angle Bhabha process, which agree 
within ~ 0.1% for integrated cross sections and 
~ 1% (or better) for differential distributions [12j 

eng. 

NNLO QED calculations are essential to estab- 
lish the theoretical accuracy of existing genera- 
tors and, if necessary, to improve it below 0.1%. 
In particular, the one-loop corrections to single 
hard bremsstrahlung should be calculated for full 
Bhabha scattering, to get a better control of the 
theoretical precision. 

For next generation e + e~ accelerators 
(ILC/GigaZ), if a 10 -4 accuracy is assumed, 
present MC Bhabha programs need to be im- 
proved by the inclusion of weak and two-loop 
QED corrections, as well as beamstrahlung and 
new, more precise Aaj^j parameterizations. 
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